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Abstract. We consider ensembles of random matrices, known as biorthogonal ensembles, 
whose eigenvalue probability density function can be written as a product of two determi- 
nants. These systems are closely related to multiple orthogonal functions. It is known that the 
eigenvalue correlation functions of such ensembles can be written as a determinant of a kernel 
function. We show that the kernel is itself an average of a single ratio of characteristic polyno- 
mials. In the same vein, we prove that the type I multiple polynomials can be expressed as an 
average of the inverse of a characteristic polynomial. We finally introduce a new biorthogonal 
matrix ensemble, namely the chiral unitary perturbed by a source term. 
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1. Introduction 

Suppose that we have a set of ./V real random variables {x±, . . . ,xn} such that their probability 
density function (p.d.f.) is given by 

Pn(xi, ...,x N ) = — det h(xj)] ij=hiN det Mxj\ j=h ^ N , (1) 

where Zn is the normalization constant. We require all variables to lie on the same interval 
JCl. Consider the n-point correlation functions 

N\ 1 f f 

Pn,N{xi, ...,x n ):= — - n y ~^ J dx n+i ■J dx N p N (xi, . . .,x n ,x n+1 . . .,x N ). (2) 

/ I 

Now assume that the matrix g, with elements gij := Jj dxrji(x)£j(x), is not singular. Then, one 
can show that the n-point correlation functions can be written as the determinant of an n x n 
matrix: 

p n ,N(xi, ...,x n )= det [K N (x i ,x j )] i j=1 n , (3) 
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where the function if at, called the kernel, is given by 

N N 
K N (x,y) = ^2 'M- r i r '\;s : . ; !.</)• ^2di,kCj,k ■= Si J. (4) 

i,j=l k=l 

Subject to a minor technical constraint on g (see 0), it is possible to construct functions 
d £ Span(£i, . . . , £at) and ifij £ Span(r/i, . . . , rj^) which are biorthogonal; that is, 

J dxipi(x)(j(x) = 6ij (5) 
l 

As a consequence, we can put the kernel in a single sum form: Kj^(x,y) = Y2i=i ^ii^div)- 
Borodin ,9, has used the expression "biorthogonal ensembles" for describing systems whose 

p.d.f. can be written as in Eq. (^Q). They have been first studied by Muttalib and Frahm 

in relation to the quantum transport theory of disordered wires 0. They can also be considered 

as determinantal point processes |24| . 

Random Matrix Theory |161 1201 125] provides many instances of such biorthogonal structures. 

First, choose rjj = a^ -1 and £j(x) = e - ^^^ -1 . Then pn(x\, . . . , xn) corresponds to the 

eigenvalue density in a unitary invariant ensemble of N x N complex Hermitian matrices, defined 

through the p.d.f. 

P JV (X)cxe- w ( x ). (6) 

In that case, the simplest one, the system is described by orthonormal polynomials {pi : i = 
0, . . . , N — 1} with respect to the weight w(x) = e~ v( - x * > ; explicitly, ipi(x) = pi-i(x) and Q(x) = 
w{x)pi-\{x). 

Second, suppose that the we break the unitary invariance of © by an external source [TT1IT2]. 

Pjv(X|A) oc e- w ( x ) +TrAX (7) 

where A is a non-random N x N Hermitian matrix. When the eigenvalues {a±, . . . , oat} of A 
are all distinct, it is possible to show that the p.d.f. for the eigenvalue of X is of the form (^Q), 
with rji(x) = x l ~ l and £j(x) = e - v ( x )+ a i x . it has been proved by Zinn-Justin |341 13 5 j that the 
n-point correlation functions comply with © . In the general case where some of the parameters 
are equal, Bleher and Kuijlaars [2] have shown that the models defined by © naturally lead to 
multiple orthogonal polynomials (see for instance [TJ 32 ). In particular, they have proved that 
the (monic) multiple polynomial of type II and having degree N, P{x) say, is simply described 
as the expectation value of the characteristic polynomial, 

P{x) = (det(zl - X)), (8) 

where the average is taken with respect to the p.d.f. (J2J). 

In this paper, we show that the multiple polynomials of type I, here denoted by Q(x), can 
also be seen as averages over perturbed matrix ensembles, 

Q{x) = Res/det(zl-X) _1 \ for z G C \ R. (9) 

z=x \ I 

In the previous equation, the residue is defined through 

fix) =: Res f dt^-. (10) 

z=x J Z — t 
I 
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We also obtain a similar expression for the kernel, 

„ r x 1 „ /det(sl-X)\ , s 

K N (x,y) = Res ( ) . 11 

x — y z=y \ det(zl — X) / 

This expression was first proposed in |22| for Gaussian ensembles (i.e., V(X) = X 2 ), and for 
general unitary invariant potentials © in |1U| . 

We also give a matrix model that possesses a new biorthogonal structure: the perturbed chiral 
Gaussian unitary ensemble (chGUE). The chiral or Laguerre ensemble plays a fundamental role 
in the low energy limit of QDC [SSI- It also appears in multivariate statistics; more specifically, 
a chiral ensemble is equivalent to the matrix variate Wishart distribution. The presence of a 
source term in the matrix model describes a non-null sample covariance matrix |25| . For more 
information on this relation, see [Hll7j. For the parameter value a = 1/2 the perturbed chGUE 
gives the p.d.f. for non-intersecting Brownian paths near a wall |26j . and similarly for a a non- 
negative integer it corresponds to non-colliding systems of 2(q + l)-dimensional squared Bessel 
processes [27.'. The biorthogonal functions of the perturbed chGUE are related to the modified 
Bessel functions of the first kind. In a special case, these functions previously appeared in papers 
by Coussement and Van Assche |131 114) . 



2. Kernel and ratio of characteristic polynomials 

For any ensemble composed of matrices having real eigenvalues, it is well known (see e.g. J 
the correlations functions can be generated by averaging ratios of characteristic polynomials: 



Pu,n{x\, . . . , x n ) = Res ... Res 

Z \ = x i z n =x n 



This formula can be proved by using 



d n /-A-det(ifcl-X) 



n 



dyi"- dy n \fj^ det(^l - X) 



(12) 

Vi=Zi 



det(yl - X)- 1 ^ det(yl - X) = tr— 1— , (13) 
dy yl - X 

and by expressing the matrix average as an average over the eigenvalues. 

In the physics literature, the residue operation is often replaced by the use of Green's functions 
and density operators. This can be understood as follows. The n-point correlation function is 
the expectation value of a product of n density operators; that is, 

Pn,Jv(»i, ...,x n ) = (p(xx) ■ ■ ■ p(x n )) , (14) 

where 

p := tr5(xl - X) (15) 

and it is assumed that the points are not coincident. For our purposes, the Dirac delta function 
has to be defined as 

5(x) = - 5 6 = -Im— ^— , e^0+. (16) 

7T X z + e Z 7T x — ie 

But the advanced Green function is given by 

6 "W :=tI (x- ie )i-x - < 17 > 

Hence 

p(x) = -ImG-(x). (18) 

7T 
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It is worth mentioning that the previous formalism allows us to rewrite Eq. (jl()j) as 



Res / dt 



fit) _ 1 

Z — t IT 



Im / dt 



f(t) 



x — ie — t ' 



(19) 



The imaginary part can be taken in two ways: 1) forming the Dirac function 5(x) inside the 
integrand, then integrating; 2) deforming the contour of integration, in order to remove the 
imaginary part from the integrand, then subtracting the imaginary part of the whole integral, 
which is equivalent to closing the contour around the single pole x. 

Proposition 1. Consider a matrix model with an eigenvalue p.d.f. given by Eq. (^Q). Let rji(x) = 
x t ~ 1 and £i(x) such that g^j := fj dxr)i(x)£j(x) defines a non-singular matrix for all i,j = 
l,...N. Then 



1 „ /det(xl-X) 
Res 



x — y z=y \ det(zl — X) / ' 
where z is a complex number with a non-null imaginary part. 



(20) 



Proof. We want to prove that the previous expression for K^j{x^y) is equivalent to Eq. (J3J). Let 
us denote the r.h.s. of (|20|) as (x — y)~ l ZZ Lj<[{x,y). From Q we have 



L N (x, y) = Res / dx\-- J dxN det [&(a?j)]^- =1 det 
I I 

By symmetry of the integrand, this can be simplified 

Ln(x, y) = N\ Res / dx\ £i(xi) • • • / dx^ £jv(x/v) det 

z =y J J 



JV 



JV 

n 



ij=l" Z—Xi 



J-l 



N 



JV 

n 



i,j=l f J -Z-Xi 
1=1 



But one proves by induction that 

JV 



JV 



JV 



f[ 1 =y - f[ 1 , 

j=1 Z \Z Xi) ^ Xj 



(21) 



(22) 



Moreover, det 



J-i 



JV 



A(xi, . . . ,x N ) = Yl 



l<i<j<N\^J 



(xj — Xi) is the Vandermonde determi- 



X j X -j 



det 



JV 



(-l) N - i A ( - l \xi,...,x N ), 



(23) 



nant. From this we deduce 

n 

where AW(a?i, . . . , xjy) = A(xi, . . . , Xj+i, . . . , xn). By substituting formulae (|2*2*)) and ((2*51 
into (|2*T|). we get 



v 



L N (x,y) = N\ V^l^Res / cfciM^l 



i=l 



dx N Cjv(sjv) 



2 X 7 ; 



V 



(x — Xj ) det 



j'=i 



j=l,...,N 
fc=l,...,JV-l 
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The two factors on the last line can be replaced by the Vandermonde 

A(xi, . . . ,Xi_i,x i+ i, . . . ,x N ,x) 

Taking the residue then gives 



N 



N 



L N (x,y) = (x-y)Nl ^(-l)^(y) [ [] / dxj^Xj) det 

<3=1 T 



1=1 



k=l,...,N 

[aJ fe-1 ]fc=l,...,JV 



Recalling rji(x) = x % , gij = J* dxr)i(x)£j(x), and integrating the determinant row by row, we 
find 



N 



Ln(x,v) = (x-y)Nl Y]Ci{y) det [g jyk ] j=i,...,n [Vj(. x )) i=1 N [9j,k] j=i,...,N 

' ■ J fe=t+l N. 



i=l 



fc=l i— 1 



We now return to K^{x^y) = (x — y) l Z N l L^{x,y). From the p.d.f. (JJ), we have 



This leads to 



Zjv = 7V!detg, g^c" 1 . 

N 

K N {x,y) = Y J Uy) detK«, 



(24) 



i=i 



where 







-(i-l)x(i-l) 
(iV-i+l)x(i-l) 



[EftLl%0)c fc 1 



-l)x(JV-i) 
iJj=i,...,AT 1 (Af-i+l)x(Af-i) 



Therefore K N (x,y) = ]Cffc=i (z)cfc,i, as desired. 



□ 



Before going to the next section, let us show that Eq. (|2U|) readily furnishes the Christoffel- 
Darboux formula for orthogonal polynomials. We define the orthogonal ensemble by r}i(x) = x l ~ l 
and £i(x) = x t ~ 1 w(x), where w(x) > is the unnormalized weight function. Let p%(x) = 
x l + cix' l ~ l + . . . denote the monic orthogonal polynomial with 

/ dxw(x)p n (x)p m (x). 

1 

We find det[a^~ ] = det[pj_i(xj)] and = N\ h n . On the one hand, by proceeding as 

in the proof Proposition 1, we find 



N 



K N (x,y)= 1 ResV] / dxi w(xi)p (xi) . . . [ dx N w(x N )p N -i(x N -i) 

x -yUi=o hn z =yfri J T J 



x - — — det 

Z Xi 



[pj-l(x k )] j=l,...,N [pj-l(x)], =1 N [pj-i{x k )] j=l,...,N 
k=l,...i-l k=i+l,...N. 



• (25) 



By virtue of the orthogonality of the p n 's, the latter equation is equivalent to 

N-l 1 

K N (x, y) = w(y) ^ T-p n {y)p n {x) . (26) 



n=0 
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On the other hand, making use of 

det Mx i )]^ j=1 = Y[w(x i )det 



N 



we get 



i=i 



Po(xi) ... PN-2(xi) 
p (x 2 ) ■■■ PN -2(^2) 



Po(x N ) 



p N - 2 {x N ) 



Z — Xl 



Z-X2 



z—x jv 



(27) 



(x 



/N-l \ 

V) [ IJ K I K N (x,y) 
\n=o J 



z=y 



Res / dxi 



Z — Xl 



(ix 2 w(x 2 )po(^2) • • • 

Po{xi) 



cIxnw(xn)pn-2(xn) det 



Po(^at) 
Po{x) 



We finally integrate the determinant row by row and arrive at 

ho ... 



(x - y)\ Y[ K) K N {x, y) = w(y) det 



,n=0 











hi 



1 pi(y) 

1 pi(x) 

The Christoffel-Darboux is established by comparing the last expression with Eq. 11261): 












Pn{x\) 

Pn(x n ) 
Pn{x) 






(28) 



h N - 2 

PN-2{y) pn-x{v) pn{v) 

PN-i{x) Pn-i(x) Pn{x) 



N 1 

Y] 7 Pn-l{y)Pn-l{x) 

^ hn-1 



1 p N {x)p N -i(y) - p N ^i{x)p N (y) 



n=l 



fov-i x - y 

3. Multiple polynomials of type I and II 



(29) 



Let us first review some properties of multiple polynomials (e.g., see ^ 132.)- These 
mathematical objects are associated to D distinct weight functions W{{x). (We limit ourself 
to the so-called AT systems in which the support I is the same for all weights i/jj.) They are 
indexed by n = (ni, . . . , no), a composition (or a multi-index) of length D and of weight N; 



that is, an ordered sequence of D non-negative integers rn such that \n\ := Y2i=i n 



N. 



i(0 



To each composition n, we associate -D multiple polynomials of type I, here denoted by A f 
where i varies from 1 to D. They are generated by the multiple function 

D 

tW^ ( 30 ) 

i=i 

This function satisfies a simple orthogonality condition 



Q H (x) = w(x) ■ A H (x) :=Y,^\x)Af{x). 



dx x 3 Qfi{x) 



0, j = 0,...,\n\ -2, 

1, 3 = H ~ !• 



(31) 



A NOTE ON BIORTHOGONAL ENSEMBLES 7 

Note that the degree of A^S is assumed to be rii — 1 (technically speaking, we only work with 
perfect systems). 

The multiple polynomial of type II characterized by the composition n is denoted by P^. It 
is a monic polynomial of weight |n| that complies with D orthogonality relations, 

dxw {i \x)x j P H {x) = 0, j = 0, ...,m-l, i = l,...,D. (32) 



l 



The multiple functions and P^ provide a biorthogonal system. Indeed, first fix ft 
(rii, . . . ,nn) with |n| = iV. Second, choose a sequence of compositions such that 



\n y >\ = i and 7V- < n- 



for all i = 0, . . . , N — 1 and j = 1, . . . , D. For instance, one could take 

n(°) =(0,0,0,...), 
n« =(1,0,0,...), 

rj(m) = (ni,0,0, ...), 
n( ni+1 ) = (m,l,0,...), 

n ,)V) = (ni,n 2 , . . . ,n/j,). 

Third, define 

Pi'-=Pnti), Qi--=Qna+v> i = o,...,N-i. (33) 

Then, we see from relations (j31|) and Q32|) that these functions are biorthogonal: 

J dxPi(x)Qj(x) = Sij for all i,j = 0, . . . , N - 1. (34) 

Proposition 2. Suppose that we have a matrix ensemble with eigenvalue p.d.f. of the form (JXJ) 
with rji{x) = or equivalently a monic polynomial of degree i — 1, and 

[tub, ■ ■ ■ ,Zn] (x) = [w^ (x),xwW (x), . . . , x™ 1 - V 1} (x), w {2) (x),xvjW ( x ), 

. . . , x n *- l wM (x),... (x),xw (D ^ (x),..., x nD ~ l wW ( x )] . (35) 

Suppose moreover that g := [gi,j\ij = \ is non-singular, where gij = fjdxrn(x)£j(x). Let 

Q H {x) = Res/det(zl -X) _1 \, (36) 

z=x \ I 

for zeC\l, and 

Pa{x) = (det(xl - X)). (37) 
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Then 



Qn{x) 



Nl 

Zff 



9i,i 



91,2 



Pti(x) = ^~ 



Zn 



9n -1,1 9n-i,2 ■ ■ ■ 

Zi(x) ... 

91,1 ■ ■ ■ 9i,n 

92,1 ■ ■ ■ g2,N 



9l,N 



9n-i,n 



(38) 



r/i(x) 



9N+1.1 ••• 9n+i,n Vn+i{x) 



(39) 



where Zn = Nl det[gij]fj =1 . Furthermore, and P^ are the only functions satisfying Eqs. 
(|3TJ) and (|32). 

Proof. In 6 , Bleher and Kuijlaars have shown that ()32|) holds true if the type II polynomials are 
given by (|37|) . or equivalently by (|39|). with w^\x) = w{x)x d ^~ l e aiX where dj = j — Y^k=i n ki 
for i such that Ylk~=i n k < j < J2k=i n k- The generalization of their argument to our case is 
immediate. So, let us concentrate on type I functions, defined in Eq. I|36|) . 

Firstly, by following the method exposed in the proof of Proposition 1, we find 

v 



Qn(x) = ^-Res j dxx--- J dx N detMx J )$ j=1 det[xy% =1 f[—}- 

T T 1=1 



Nl 

— Res / dxiii{xi) 

Z]\[ z=x 



dxN £n (xn) det 



[r]i(xj)]i=i...N-i 
j=l,...,N 

MM 

lz-Xj J j=l,...,N 



Nl f f 

Res / dx\ ■ ■ ■ / dxjydet 

Zn z=x J J 
i i 



\x) 1 ^j(x)] i=1 ... N -i 
:i i v 



2 — Xj }j = l....,N 



The last line obviously leads to Eq. (JSHJ). 
Secondly, we choose rn(x) = x l ~ l and set 



dx Wi(x)r]j(x), 



(40) 



so that the determinantal expression of becomes 



Qfi(x) 



Nl 

Zn 



i=l,...JV-l 



I'H+j-lh 



\i=l,...,N-l 
j=l...,ni j=l...,n,2 

[rjj (x)w x (x)] j=l ^ ni [Vj (x)w2 (x)] j=h 



n 2 



Hence, 



, Nl 

dxx Qn{x) = — 

ZjN 



[^i+j_l]t=l,...,JV-l 
j=l...,m 



+ j-l\i=l,...,N-l 
j=l...,n 2 



[Vj(x)w D (x)]. =h ni 



i=l,...,N-l 
j=l...,n D 



(41) 



+j-lH=l,...,N-l 
j=l...,n D 



[h™ 



'j+klj=l,...,n D 



j+klj=l,...,m L j+kl j=l,...,n 2 

We see that the r.h.s. is null when the last row equals one of other the rows; i.e., when k = 
0, . . . , N — 2. For k = N — 1, the determinant simply becomes det[5ij]^- =1 = Zjq/Nl. This 
completes the proof of the orthogonality condition (|31|) . 



A NOTE ON BIORTHOGONAL ENSEMBLES 



9 



We finally show the uniqueness of the type II multiple function. Expression (|38j) tells us that 



Qn(x) = E 



i=l 



(0. 



where ^4 



(i) 



(O^rij — 1 _|_ c w x m 



~ 2 + . . . + Cn}- This means that, 



in order to determine uniquely, we have to fix the N = \ft\ coefficients cjp (j = 1, ... 
i = 1, . . . ,D). But Eq. (|31|) furnishes exactly N linear equations. The matrix for that linear 
system is 



g 



9l,i 



9N,1 



9l,N 



9n,n 



+ j_lJi=l,...,JV 
j=i...,m 



+ j-l\i=l,...,N 
j=l...,n 2 



14 D) 



+j-l\i=l,...,N 
j=l...,n D 



But by hypothesis detg ^ 0. Consequently, the solution for the coefficients, and therefore Q^, 
is unique. □ 



4. Perturbation of chiral unitary ensembles 

As we mentioned in the introduction, non-trivial biorthogonal matrix ensembles exist. For 
instance, ensembles of N x N Hermitian matrix with a p.d.f. proportional to exp(— trV(X) + 
trAX), where A is a fixed N x N Hermitan matrix, naturally lead to biorthogonal systems. 
Specifically, let £i(x) = exp(— V(x) + a^x) and let a := {ai, . . . ,ajy} denote the eigenvalues of 
A. Suppose that some of the Oj's coincide, i.e., 

a ni —> Oni-l ai = 6i, a ni+ri2 — > a n2+ni -i a ni+ \ = b^, (43) 

and so on. Symbolically, this is written as 

a = b™. (44) 
Then from the definition of the function £j one sees 



n 

n>0 



uniformly for \dj — Oj| < 00, so that [§] 

im rw,<7v(«;-ai) uliUT^miKkKiKDibi-h)^ 

where [£i,...,£jv] is given by Eq. ()35[). with Wi(x) = exp(— V(x) +bix). Therefore, these 
biorthogonal ensembles can be studied with the help of multiple polynomials (see Proposition 
2). In this section, we provide other examples of biorthogonal matrix ensembles. This time, 
weight functions are of the form w <yt \x) = x a e- v ^I a {2^fxbl) 1 where I a is the modified Bessel 
function (see below). 

Suppose a := M — N > 0. Let X = pQj], be a random M x N (non-Hermitian) complex 
matrix drawn with probability 

P(X)(dX) = e - trV/(xtx) e Rc(trXAt) (dX). (47) 
A is a fixed M x N complex matrix and (dX) denotes the normalized volume element of C MxN , 



M N 

(dX) := ^Hl[dRe(X itj )dIm(X h3l 

i=ij=i 



Xr 



U e U(M), V e U{N). (49) 
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The "potential" V has to be chosen in a way that guarantees the positivity of trV(X). When 
A = OmxN and V(x) = x, the p.d.f. (jITjl defines the chiral Unitary Ensemble (chUE), which is 
simply related to the Laguerre Unitary Ensemble (LUE) (for more details, see |20| Chapter 2]). 

We now want to get the eigenvalue p.d.f. associated to (|47jl. This can be done through a 
singular value decomposition of X: 

X = UXrjVt, (48) 

where 

diagfsi, . . . ,s N ] 

Note that the singular values s±, . . . , sjv are real and non-negative; they are the positive square 
roots of the eigenvalues of the N x N matrix X^X. A similar decomposition is possible for the 
non-random matrix; that is, A = UAdV^ with A^ = [diag[ii, . . . , ijy] Onxo] and t±, . . . , i/\r > 0. 
For the moment, we assume ti ^ tj for i ^ j. For notational convenience we set 

sf = x h t-=4oj, i = 1, . . . , JV. (50) 

By considering the transformation (|48|). and its Jacobian, we get an integral form for the eigen- 
value p.d.f. 

N 

p N (xi,...,x N ) oc JJxfe" y(:ri) Yl {xj-Xif 

i=l l<i<j<N 



(VUV) J (Ut(iU)exp jRettrXoVtApjU)}, (51) 



x 

U(N) U(M) 

where (U+dU) is, up to a multiplicative factor, the Haar measure on the unitary group U(M) 
(and similarly for V). The integration can be realized by making use of a Itzykson-Zuber type 
formula [23 EU: 



(VUV) J (UtdU)exp|Re(trX D V t A{ ) U)} 



U(N) U{M) 

c fr 1 det[j Q (2^x7)]^ =1 

' fj[ (a,iXi) a / 2 Ajv(ai, . . .,a N )A N (xi, . . .,x N ) : 

where Cm,n is a constant independent of the Xj's and a^s. Recall that the modified Bessel 
function of the first kind is specified by 



c 



{()} 



where it is assumed that a£Z, and where Ctq} stands for a positive contour that encircles the 
origin. It can be expressed as a hyper geometric function as well, 

z a/2 

I a {2z 1 / 2 ) = ——— F 1 (a + l,z), aeC, |arg(z)|<vr. (54) 
T{a + 1) 

Combining the few last equations, we obtain the eigenvalue (or singular value) p.d.f. 



, , 1 tt a -v(xi) TT fxj-xA \ Fi(a + l,aiX^ l X 
p N {x x ,...,x N ) = — r-jjafe V(Xl) [] [— ) det 



N i=l l<i<j<N ^ J 1 



T(a + 1) 



(55) 



A NOTE ON BIORTHOGONAL ENSEMBLES 



11 



L'Hospital's rule provides the appropriate density when some of the aj's coincide. Clearly, Eq. 
1)55(1 is of the biorthogonal form, with 



rji(x) = x l + lower terms, £,i(x) 



X a e -V(x) 



\Fi(a + l,aix). 



(56) 



r(a + i) 

As a consequence, the correlation functions satisfy p n ,N{xi, ■ ■ ■ ,x n ) = det [KN(xi,Xj)]™j =l and 
the kernel is given by (x — y)Kj^(x, y) = Res (det(xl — X)det(zl — X) -1 ). 

z=y 

When we perturb ensembles of Hermitian matrices by a source term, the multiple functions 
and Pfi can be obtained through Proposition 2. In that case, the composition n = (n\, . . . ,no) 
gives the multiplicity of the eigenvalues ,bo) (see limit (031)) • One might be tempted 

to conclude that this relation remains the same in perturbed ensembles of rectangular complex 
matrices. It is true that Proposition 2 still holds. However, the link between the multi-index n 
and the eigenvalues (b\, ...,&£>), or equivalently between the function £j and the weight functions 
Wi, is more involved. The following lemma and proposition aim to clarify the situation. 

Lemma 3. Let £i(x) = w a (x,ai), where 



X a e —V(x) 

w a (x,ai) := — — ——oF 1 (a + l,dix). 



fall 



Consider the limit (|43|) . Then, 



lim 



det 



it; 



(i> 



(2) 



Ui^jKN^j - oi) nf=i n;^ 1 $ rii<*</<y>(fc - &*) 



n fe ra £ 



(57) 



(58) 



where 



w a {xi,b k ) w a+ i(xi,b k ) 



W a+nk -l(Xi,bk) 



=1,...JV 



Proof. First, we suppose that, as a n — ► a n _i 



ai = 61 , the following equation holds true 



G n := lim JJ (aj - ai) 1 det[^ i (x j )}^ j=1 = 
l<i<j<N 
n-1 N 

Y[(k\r l H (a,, - bxT n J] (aj - ai)' 1 x 

k=l i=n+l n+l<i<j<N L ~- J ' " ,J, ="+ 1 .-." i=l,...,N 

Second, we consider the series expansion of I a , given by Eq. (|53[). from which we deduce 

(Oj - Oj) 



[Wa+i-tix^bi)].^^ 

[w a (xj,ai)]. = 



w a (x,aj) = ■ 

uniformly for |aj — aj| < 00. We thus have 
lim TT (aj - ai) _1 det[^(xj)]^ =1 

0-2 — >ffll 



-w a+ ^(x,ai) 



• (59) 



(60) 



l<i<j<N 



N 

n(a» - £>i)~ 2 TT (aj - a;) -1 lim (a 2 - ai) -1 
i=3 3<i<j<N 

[w a (xj,ai)] 

[w a (xi,ax) + (02 - ai)ui Q+ i(xi,ai) + O ((a 2 - ai) 2 ) ] 
t' u;Q! ( a: i' »)]i=3,...,JV 



j=l,...,Af 
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But we can subtract the first row from second without affecting the determinant, so that 



G 2 



lim 



a 2 ^a 1 =b 1 ril<i<j<Ar( a j _ a i) 



N 



ri(« 2 



n («, 



a,- X 



i=2 



3<i<j<N 



W a (xj,bi)} 

w a+1 (xj,bi)] 

(Xj,0i)] 



i=3...N 



j=l,...,N 



This shows Eq. ()59|) for n = 2. The general n case is established by induction: we return to 
(|59j): we apply (|60|) once again, i.e., 

n-l N 

lim G n = n^O" 1 II (a n -a t )- n - 1 \\ i^-^y 1 

Q'TI — I— 1 

fc=l i=n+2 n+2<i<j<N 



lim (a n+ i - 6i) 



[^0+4-1(^,61)]^^ 

[ w a(xj,Oi)] i=n+2 N 



j=l...N 



-Row(n); 



we manipulate the rows as, 

Row(n + 1) -> Row(n + 1) - Row(l) - (a n+1 - MRow(2) ... 

(n — 1)! 

and we finally get 

n AT 

n^r 1 n n 

fe=l j=n+2 n+2<j<j'<A r L " J ' "' J j=l...N 

which is G n +i, as expected. The general formula (|58|) is obtained by taking D limits similar to 



(Oj - aj J x 



[ W «+*-l( X i' 6 l)]<=l,...,n+l 

[u; a (ay,ai)] i=n+2 ,.. .,n 



Proposition 4. Consider the functions £j and w a as defined in the previous lemma. Let w a (x) 
w a (x,0). Suppose 



Then, 



ai = 6 > 0, 



a r+ i = 0, 



(61) 



lim 



detfc(xj)]& =1 

ni<i<j<Af( a j _ a «) 



nu j (-&)-nSV* 



n£ 1 1 «n£ r - 1 - 



i=l,...,jV-r 



x) w a+ i(xj,b) 



(62) 



j=l,...,JV 



Proof. We have from the previous lemma and from w a +i(x,0) = x l w a (x)th&t 

x 1 f l w a {x j ) 



lim 



dfitfcfa)]j^ 



i=l,...,7V-r 
i=l,...,JV 



to a +i_i(a;j,6) 



i=l,...r 
j=l,...,JV 
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Now, for a > — 1, it is known that 

x -Fi(a + 3,x) = (a + l)(a + 2)( i ? i(a + l,x) - F 1 (a + 2,x)\ 
This implies for b ^ and k EN, 

bw a+k (x,b) = xw a+k - 2 -{a + k- l)w a+k -i(x,b). 
The latter identity allows us to write 

w a+k (x) = 2 w a+2( ^k_yk^{x) + c a {x;b,k) 

where c a is a linear combination of w a and u> Q +i with coefficient depending on b and fc. The 
desired result is obtained by using the latter equation and by exploiting the antisymmetry of 
the determinant under the permutation of the rows as well as the invariance of the determinant 
under the transformation Row(i) — > Row(i) + Ylj^i c i,j^ ow U)- ^ 

The last proposition implies that each limit of the form a ni — * a ni -i ■ ■ ■ — > a i = b > 
give rise to two functionally independent weight functions, i.e., w a {x,b) and w a+ \{x,b). When 
we have D similar limits, with b\ > 62 > • • • > say, we get 2D weight functions if bp > 0, 
and 2D — 1 weight functions if bp = 0. 



5. Chiral Gaussian Unitary Ensemble with a source term 
In the next paragraphs, we focus on the perturbation of the chGUE. So we choose 



V(x)=x, Vk (x) = (-l) k -\k-l)!L^ 1 (x), U*) = Xae ~ X f ( l{a ^ aiX) (63) 

I (a + 1) 

where denotes the (associated) Laguerre polynomial of degree k. 

Proposition 5. Let a = {a\, . . . ,ajy} E (0, 00)-^. Then the kernel of the perturbed chGUE, as 
defined by Eqs. (|47|) and ([63)1 . is given by 

V) = £^ [ du „« e - F l{ a + 1, -xu) [ + l,-yv) jj u + a, 

r M l ~ l 

°— a 

where C_ a denotes a counterclockwise contour encircling the points —a%, . . . , — a^v frwi 
point u. Equivalently, 



W a (x) „ , , 

K N {x,y) 



w a (y) 



to I *L I *0L I *L e ^ e "" /w+ ^ e ^ fir frlLtg, (65) 

' 27ri 7 27ri 7 27ri u — v zw \v/ \ws . r; v + a r - 



C M c {0} c {0} 

b 



i=i 



where C| j is a positive contour around the origin. When some of the parameters ai 's are null, 

du is understood as lim e _> + f° 



the previous equations remain valid if f^°du is understood as lim e ^ + / e °° du. 
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Proof. In the first steps, we assume that all the elements of b are distinct. As mentioned in the in- 
troduction (see also Proposition 1), the kernel can be written as Kn(x, y) = J2^j=i ^( x ) c i,j£j(y) 
with [c jti ] = [gi,jY l (c* = g" 1 ). Explicitly, 

N 

y a e~y 



K N (x,y) = T(a + 1) E " l)\L?_ xCiij0 F x {cL + l,a jV ) 

N 00 

= ^ + 1)2 E(" 1 ) i ~ lc M J duu a+i - 1 e- u F 1 (a + l,-xu) F 1 (a + l,a j y) (66) 



where we have made use of the formula (cf. |31l Eq. (5.4.1)] and Eq. (|54|) above 



00 

L°(x)= f [ duu a+n e~ u F 1 (a + 1,-xu). 
n\T(a + 1) J 



(67) 



We now aim to eliminate the coefficient Cjj. In our case, 



00 

9i,j = ( ~r( a |\^ 1)! / dxx a e- x L?_i(x) oF^a + l,a jX ). (68) 


This can be evaluated exactly: 

.'A.., ") (69) 
As a consequence, Cjj must comply with e ak Y^iLifakY* 1 °ij = $j,k- Thus 

N N 

By substituting the last equation into (|66|) . we find 

a — y+a: ^ ^ AT t 

K N (x,y) = * \ / duu a e- u oF 1 (a + l,-xu)Y,oF 1 (a + l,a j y)e-^ll^-^. (71) 
1 (a + i) J aj ci£ 

1 m 

Clearly, the sum can be written as the addition of all the residues at v = ai (£ = 1, . . . , N) of 

N 1 -, 
^ / n rr u + dp 1 
oFx(a + 1, vy)e- v TT . 

±J - V — CLf U + V 

1=1 - 

This proves Eq. (|64j) when all a/s are distinct. The validity of this expression for the general 
case is established by continuity. The quadruple integral representation of the kernel is readily 
obtained by comparing Eqs. I|5HJ) and (|54j) . □ 

Multiple functions associated to the modified Bessel function of the first kind have been re- 
cently introduced by Coussement and Van Assche 13] . Their work, which was not motivated by 
Random Matrix Theory, corresponds to our D = 1 case. Here we provide integral representations 
of the type I and II multiple functions for D > 1. 
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To avoid any confusion between the number of weights and the number of distinct eigenvalues 
in the perturbation matrix, we have to introduce some new notations. Suppose that a = b m , 
where m = (mi, . . . ,m d ), \rh\ = N, and 61 > 62 > • • • > b d > 0. Set n = {n\, . . . ,njj) with 

mi + 1 



^3 



2 

m2 + 1 



mi 
m 2 



UD-l 



m d + l 



n D 



m d 
L 2 



(72) 



where D = 2d. If b d = 0, it is understood that = m d and D = 2d — 1. The multi- index n 
gives the correct multiplicities of the weight functions (see Section 3); that is, if bd > 0, 

w(x) = w a (x,bi), w a+1 (x,bi), w a (x,b 2 ), w a+ i(x,b 2 ), ■ w a (x,b d ), w a+1 (x,b d ) (73) 

or, if b d = 0, 



wlx) 



w a (x, bi), w a+ i(x, 61), w a (x, 62), w a+ i(a;, 6 2 ), . . . , io a (:c, 6 d ) 
where w a stands for the function defined in Eq. (|57|1 with V(x) = x . 

Proposition 6. Following the above notation, we have that the multiple function of type I is 

dv e~" oF\{a + 1, xv) 



(74) 



Qn(x) = w a (x) 



2^ n£i(«-oi) 



T(a + l)w a (x] 



dv r dz z <x-i e v(x*-i)+y* 
, ™J ™ Uli(v- ai ) 



(75) 
(76) 



while the multiple polynomial of type II is given by 



r(a + i)7 



T(a + 

(-1) 



f dw e -<^+i)+i/w ^ 

, du I \ \[u + 

T(a + l)w a {x) J J 2vri fl v 
c {0} l - x 



(77) 



(78) 



Proof. We start with the determinant representation of the type I function given in Proposition 
2. Here, 



A/ 



Ztv = Z' N A( ai , ...,a N ) = N\ ~[\e a >A( ai , a N ). 



Eq. (|3*5|l then implies 



Qn(z) = X) ff 



i=l 



V &(*)e" a ' . . A e" * 0^1(0 + 1,0^) 



W a {X) > pp- 

which turns out to be equivalent to the proposed formulae. 
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The proof for the type II function is a bit more tricky. Firstly, Eq. (j39() gives 



Pn(x) 



("I) 



N N 



A(oi, ... ,a N 



n=l 



J-l 



J-l 



i=l,...,n— 1 
j=l,...,N 

i=n+l,...,N 
j=l,...,N 



Secondly, we introduce the elementary symmetric functions, denoted by e n , and defined via the 
following generating function: 

N N 

Y[(t + <h) =: ^2 t n e N - n {ai, ■ ■ cln)- 



i=l 



Then a few manipulations give 

1 



A(oi, . . . ,ojv) 



n=0 



i=l,...,n— 1 
j=l,...,N 

i=n+l....,N 
j=l,...,N 



ejy+i-n(ai, • • • ; «tv) 



Hence 



v 



Pjj(a?) = (-1)^ '^2(-l) n rj n+1 (x)e N - n (a 1 , . . . , ojy). 



(79) 



n=0 



Thirdly, we exploit the integral representation (|57|) to obtain 

AT OO 

Pn(x) = L~ \ e - V / di^e^oFiCa + l.-artJ^ejv-nCai,... 



Oat, 



The first integral representation of Pjj is finally obtained by reconstructing the generating func- 
tion of the elementary symmetric functions. The last double integral expression is a mere 
consequence of Eqs. ((55)1 and ((54)1 . □ 

The above multiple functions can be considered as multi-parameter generalizations of the 
Laguerre polynomials. Indeed, from the integral representations of the Laguerre polynomials, 



L%(x) 



dw e 



2vri w N+1 



[1 + w) 



N+ 



m 



dw e a 



2vri w N+a+1 



(w-1) 



N 



C{o} 



'{0} 



when a £ Z, one can show that 



(-1)1*1-1 



lim Q n {x) = ^ /_ x a e- x Lf Hhl (x) 



a^O 



(|n| + a 



and 



\unP n (x) = (-l)W\n\\L? f . l (x). 

a— >0 1 

Note that the last identity is more easily shown by using Eq. ((79(1 . i.e., 

N 

Pft(x) = (-l) Ar ^n!eAr_ n (ai,...,a A r)L"(x), 



(80) 

(81) 
(82) 

(83) 



n=0 



which is valid for all a%, . . . , ajv- Note that the formulae ((75)1 and 1(77)1 furnish an alternative 
method to prove some of the properties found in |13j for the multiple polynomials with D = 1. 
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Corollary 7. Let 01 > a 2 > . . . > a N > 0. Let also n<°> = (0, 0, . . . , 0), = (1, 0, 0, ... , 0), 

re^ 2 ) = (1, 1, 0, ... , 0), and so on till n( N ' = (1, 1, . . . , 1). Define Pi = P-(«) and Qi = Q^i+i) ■ 
Then 

N-l 

K N (x,y)=Y J Pi(x)Qi(y)- (84) 
i=0 

Proof. We substitute the formula |17| 

U + di 1 ^ lli=l ( U + 



_1 "I r U + Oj _ 1 

?) 11 ?! -L li 91 ^- ^ 



u _ „ ii, + ai u - v ^ nLi(« + «*) 

in (|64|). We then use the fact that 

e 1 ' o-^i (a + 1, — 



2tti u — v 



and obtain 



^r(a + l) 2 7 AX _/ 2m Y[ k ( v + a 

The comparison with Eqs. (|75|) and (|77|1 finishes the proof. □ 

The general Christoffel-Darboux (CD) formula involving multiple polynomials of type I and II 
has been found by Daems and Kuijlaars JH]- The complicated relation between the weights and 
the perturbation eigenvalues makes difficult the extraction of the general CD formula from the 
integral representation of the kernel. This is in contradistinction with the perturbed Laguerre 
ensemble in which the CD formula is readily derived from integration by parts [HI El 



6. Concluding remarks 

The integral representations of both the kernel and multiple functions of type I and II provide 
tools for studying the asymptotic behavior of the perturbed chGUE. Of particular interest are 
the finite rank perturbations [UEHED]- in such systems, the eigenvalues of the perturbation 
matrix A satisfy 

a= (ai,...,a r ,0,... ,0). 

The rank of the perturbation is given by r; it is finite in the sense that limjv^oo N _1 r = 0. It 
has been observed in that the kernel of the perturbed Laguerre ensemble can be decomposed 
as a unperturbed kernel Kn_ t , where 

K N (x,y) = N] ^ (L%_ 1 {x)L%{y) - L%{x)L%_ 1 (y)) , 
{Jy + a — 1)1 x — y V / 

plus a sum of r projectors Pi® qi, where pi and % are respectively related to the type II and I 
multiple Laguerre functions. A similar decomposition exists in the perturbed chGUE; explicitly, 

r 

K N (x,y) = K N - r (x,y) + ^p fc (x) q k (y), 

k=i 
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where 



fc-i 



,"/, (•'') - — / duu N+a - r Y[(u + ai)e- u oF 1 (a + l,-.ra) 

T(a +1) J t~ 

ii « — J- 

dv e v oFi(a + 1, —xv) 



r(a + l) J 27ri v JV-rnf =1 ( u + ai ) 

k— a 

The asymptotic correlations of the chGUE with perturbation will be considered in a forthcoming 
paper. 
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